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We consider fretting wear due to tangential oscillations of two contacting bodies. For small oscillation 
amplitudes, the wear occurs only in a circular slip zone at the border of the contact area. With increasing 
number of cycles, the wear profile tends to a limiting form, in which no further wear occurs. Under 
assumption of a constant coefficient of friction, the limiting form of the wear profile does not depend on the 
particular wear criterion and can be calculated analytically. An explicit analytic solution is presented for 
arbitrary initial shape and illustrated for the cases of parabolic and conical shapes. 

Fretting wear occurs if two bodies are pressed against each other and subsequently subjected to tangential 
oscillations with small amplitude. It is of importance for any technical system with joints subjected to 
oscillations. Very much attention was paid to fretting in such applications as fretting of tubes in steam 
generators and heat exchangers 13 , joints in orthopedics 4 , electrical connectors 5 , and dovetail blade roots of gas 
turbines 6,7 as well as many others. Depending on particular properties of materials and on loading conditions, 
fretting wear can lead either to a progressive wear or to some final state in which no further wear occurs. The proof 
of the existence of the final state was done by Ciavarella and Hills 8 for the case of an arbitrary two-dimensional 
(plane strain) contact problem. In the present paper, we show that their conclusions are also valid for three- 
dimensional contacts of axis-symmetric bodies and derive the shape of the final, no-wear state. 

Consider a rotationally symmetric profile which is brought into contact with a rigid surface and then oscillates 
in tangential direction with a given amplitude uf \ Under assumption of sticking condition, the tangential stress 
has a singularity at the border of the contact area 910 . In a contact with any finite coefficient of friction, this will 
lead to the appearance of a circular slip zone at the border of the contact area 1113 which leads to wear in this 
zone 8 . In the present paper, we will assume that the friction can be described by a local formulation of the 
Amonton's law: The surfaces in contact are in the sticking state if tangential stress x is smaller than normal 
pressure p multiplied with a constant coefficient of friction ji, and the tangential stress remains constant after the 
onset of sliding: 



t < f.ip, stick 
t = fip, slip 



(1) 



At the circular border of the stick region with radius c, the critical condition t(c) = ^ip(c) is fulfilled. Inside this 
region, the condition t < up is valid. Due to wear outside of the sticking region, the local pressure in the sticking 
region will increase and outside decrease further, independently of whether the experiment is done under 
conditions of constant normal force or constant indentation depth d. This will lead to a progressive wear outside 
of the region of stick. The wear process will advance until the pressure in the sliding region becomes zero. In this 
limiting state, the inner parts of the contact will still remain in the sticking state, while the wear rate in the outer 
parts of the contact tends to zero. The final state of no wear can be considered as a sort of "shakedown" state, in 
which no further inelastic processes occur. The detailed kinetics of the profile depends on the wear criterion used 
as well as on the loading conditions (controlled force or controlled indentation). In the most cases, the Reye- 
Archard-Khrushchov wear criterion is used, stating that the wear volume is proportional to the dissipated 
energy 1416 . According to this wear criterion, the wear rate vanishes if either the relative displacement Au x of 
the bodies or tangential stress in contact is zero. In non-adhesive contacts, the latter means vanishing of the 
normal pressure p. The no-wear condition thus reads: 



No wear condition 



either p = 0 
or Am., = 0 



(2) 
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From these conditions, it follows that the pressure in the final state is 
non-zero only inside the stick area and vanishes outside. The whole 
process and the condition (2) for the final state were formulated in 
detail by Ciavarella and Hills in Ref. 8 . They suggested that it should 
be possible, using this condition, to deduce the final shape. This is 
done in the present paper. 

Results 

Analytical solution for the final state with the method of dimen- 
sionality reduction (MDR). To find the final profile we use the 
method of dimensionality reduction (MDR) 1719 . This method 
allows solving the normal and tangential contact problems for 
axisymmetric bodies by mapping them to a one-dimensional 
contact of properly defined elastic foundation. We would like to 
stress, that in spite of analyzing a contact with one-dimensional 
elastic foundation, the MDR provides exact solutions for the initial 
three-dimensional contact problem. The corresponding proofs can 
be found in Ref. 19,20 as well as in Chapters 17 and 18 of the 
monograph 18 . In the following, we shortly recapitulate the basics of 
the method of dimensionality reduction. 

The main steps of the MDR are the following. Given a three- 
dimensional profile z =f(r), we first determine the equivalent one- 
dimensional profile 20 



g(x)- 



f'(r) 



\fx 



.dr. 



The back transformation is given by the integral 

*(*) 



-.Ax. 



(3) 



(4) 



The profile (3) is pressed to a given indentation depth dinto an elastic 
foundation consisting of independent springs with spacing Ax whose 
normal and tangential stiffness is given by 



k z =E*Ax 
k x = G*Ax , 

where E* is the effective elastic modulus 

1 1 - v? 1 - v 2 

— = 1 H 

E* Ei E 2 

and G* the effective shear modulus 

1 _ (2-vQ (2-v 2 ) 
G* 4Gi 4G 2 



(6) 



(?) 



Ei and E 2 are the Young's moduli, Gi and G 2 the shear moduli of 
contacting bodies, Vi and v 2 are their Poisson-ratios. Note that 
throughout this paper, we assume that the contacting materials sat- 
isfy the condition of "elastic similarity" 

l-2vi l-2v 2 



Gi 



(8) 



guaranteeing the decoupling of the normal and tangential contact 
problems 9 . In particular, this condition is always satisfied in the 
important case of the contact between a rigid body and an incom- 
pressible elastomer (both sides of equation (8) are then zero). 
The resulting vertical displacements of springs are given by 



u z (x) = d-g(x) 

and the linear force density 

q(x) = E* u z (x) = E* (d-g(x)) . 



(9) 



(10) 



The contact radius a is given by the condition 
g(a) = d. 



(11) 



According to MDR rules, the distribution of normal pressure p in the 
initial three-dimensional problem can be calculated using the follow- 
ing integral transformation 18,20 : 



it 



\/x 2 — i 



2 n 



\/x 2 — r 



.Ax. 



(12) 



If the profile is moved tangentially by u x °\ the springs will be 
stressed both in the normal and tangential direction, and the radius 
c of the stick region will be given by the condition that the tangential 
force k x u x 0 ^ is equal to the coefficient of friction fi multiplied with the 
normal force k z u z (c): 



C?uf = nF{d-g{c)). 



(13) 



Let us denote the initial three-dimensional profile as /o(r), the 
corresponding one-dimensional image as go(x) and the limiting 
shakedown shapes as f x (r) Midg x (x) correspondingly. As discussed 
above, the pressure outside the stick area must vanish in the limiting 
shakedown state: p(r) = 0, for r > c. From (12), it follows that 



g'(x) = 0 and g{x) = const = g 0 , for c < x < a. 



(14) 



From the condition (11), it then follows that the one-dimensional 
profile in the shakedown state has the form 

go( x )> f° r 0 <x< c 



d, for c < x < a 



(15) 



This shape is schematically shown in Figure 1. The three-dimen- 
sional limiting shape can now be calculated by the back transforma- 
tion (4): 



f a {r), for 0<r<c 



(5) U(r)- 




^Lchc+^d 

x 2 n 



-. dx, for c<r <a 



(16) 



Examples of limiting profiles are shown in Figure 2. 



Limiting shapes for indenters with different initial shape. I. 

Parabolic indenter. The initial profile is in this case/o(r) = r 2 /(2R), 
and the corresponding one-dimensional MDR-image g 0 (x) = x 2 /R. 
The radius of the stick region is given by the condition (13): 



(0) 



R\ d 



G M; 

E*~li 



(17) 



According to (16), the limiting three-dimensional profile has the 
form 




Figure 1 | One-dimensional MDR-image of the final "shakedown" 
profile. 
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Figure 2 | Examples of limiting profiles with c = 0.5 for initial (a) parabolic form and (b) conical form. Dimensionless variables are defined in Eq. (19). 



— , for 0<r<c 
2R 

d ( d | arcsin (-\ \ 1 1 — ( - 1 , for c < r < a 

n V 2RJ r kR \rJ V \rJ 



an accuracy not less than 15%. Using this approximation, one can 
easily calculate the dependence of the outer wear radius a on the 
(18) inner wear radius c: 



Normalizing all vertical coordinates by the indentation depth d and horizontal 
coordinates by the contact radius of the initial profile, a 0 = VRd, 



f=f/d, d=d/d=l, 

r = r/o 0 , x = x/ao, c = c/a 0 , a = a/a 0 , 



(19) 



we can rewrite these equations in the dimensionless form 



'~r 2 r 
— , for 0 < r < c 
2 



20) 



2 ( f 2 \ e rc fcV r . 
1 1 arcsin \ / 1 — = , for c<r<a 



The non-dimensional form of the limiting profile thus depends only 
on one parameter 0 < c < 1. 

The contact radius, and thus the outer radius of the wear region, is 
given by the condition / x , (a) =fo(a): 



2 i a \ c ac 

1 1 arcsin \ j 1 ■ 

\ 2 ) a n 



c\ 2 a 2 



(21) 



In the limiting case 1 = 0, the contact radius becomes a = y2. 
The total force can be calculated as 



Fn = 2 



E*(d-g(x))dx = 2 
or, under consideration of (17), 



E*(d-x 2 /R)dx = 2E* ( dc- — \{22) 



F N =-E*R 1/2 



E* H J V 2E * V 



(23) 



Profiles (20) are shown in Figure 3 for a representative set of 
parameters. 

For c/r~l, the exact solution (20) can be expanded in series: 



/«(r>" 



for 0<?<c 



(24) 



2 | J — - J. for c<r<a 



c\ 2 c 



(25) 



The normalized wear volume in the limiting state can be evaluated 
with the same accuracy as 



V rx ~7t(c-1) 2 , 
where y x , = Vaa/idc^), and V x is the wear volume. 



(26) 



II. Conical indenter. The initial profile is in this case f 0 (r) — r tan 0, 
and the corresponding one-dimensional MDR-image go(x) = 

— \x\ tan 6. The radius of the stick region is given by the condition (13): 



n tan 0 \ E* f.i 

The limiting three-dimensional profile has the form 
( r tan 0, for 0 < r < c 



(27) 



r — \lr 2 — c 2 \ tan 9 + d [ 1 arcsin V for c < r < a ^ 



UAr)- 



Introducing the same dimensionless variables (19), but normalized by the 
2 d 

initial contact radius a 0 = for the cone, we come to the dimen- 

7i tan 0 



sionless presentation 




While this equation is asymptotically exact for c/r-»-l, it gives a good Figure 3 | 3D profiles in the final state according to equation (20). 
approximation in the whole definition range of variables i and r with Parameters: 9 linearly increasing c from 0. 1 to 0.9. 
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r 

Figure 4 | 3D profiles in the final state according to equation (29). 

Parameters: 9 linearly increasing c from 0.1 to 0.9. 



boundaries, at which the wear completely disappears and the normal 
stress becomes infinitely large, the above conditions are surely not 
applicable in the immediate vicinity of the boundary of the stick 
region. It is an interesting task for the future research to clarify the 
influence of these factors on the fretting wear and the very existence 
of the final state. 

Methods 

For analytical solution, we used the so-called Method of Dimensionality 
Reduction (MDR). The MDR is a simple technique of solving contact problems 
for axi- symmetrical bodies. In this method, a three-dimensional contact problem 
is replaced by a contact of a transformed profile with a properly defined elastic 
foundation. The most comprehensive description of the method with all necessary 
proofs is given in the recent monograph 18 and in the dissertation 19 . The basics of 
the method are further described in 17,20 . The MDR provides the exact three- 
dimensional relations for indentation depth, contact radius and normal force as 
well as radius of the stick region and the tangential displacement for the tangential 
contact. Further, it allows restoring the correct three-dimensional distributions of 
normal and tangential stresses in the contact area, thus providing the complete 
solution of the three-dimensional problem. 



r, for 0 < r < c 



— (i — Vr 2 — c 2 ^j + ( 1 arcsin ( -) ) . lor r < r < a 



(29) 



The contact radius a is given by the condition /qo (a) =fo (a), from 
which it follows that 



71 I C 

arcsin - 

2 la 



(30) 



In the limiting case c = 0, the contact radius becomes a — n/2. 
The total normal force is equal to 



F N =2 



E*(d-g(x))dx = 2 



E*\d — — x tan 0 )dx 



(31) 



--2E !cfc--c 2 tantf 



or under consideration of (27) 
IE* 



Fn = 



n tan 0 



G*zf> 
E* f.i 



(32) 



The profiles (28) are shown in Figure 4. 



Discussion 

In the present paper we considered wear of axis -symmetrical profiles 
due to tangential oscillations with small amplitude. Following former 
suggestion of Ciavarella and Hills, we have proven that there exists a 
final state in which the wear rate tends to zero. Using the method of 
dimensionality reduction, we have found an analytical solution for 
the shape in this final state in the closed integral form, Eq. (16). 
Explicit integration has been carried out for parabolic and conical 
profiles. The governing parameters of the final shape are indentation 
depth and radius of the stick region. We would like to discuss the 
main assumptions of this paper which affect the practical applicabil- 
ity of the theory: (a) using the Coulomb's law of friction with a 
constant coefficient of friction, and (b) assuming a local wear law, 
stating that the wear at each point does depend only on the pressure 
and relative tangential displacement at this point. Both assumptions 
are widely used but are known not to be completely correct: the 
Coulomb friction law is only approximately valid if the normal stress 
is not too small or too large 10 and the derivations of the Khrushchev - 
Archard wear relation show that the wear process is in reality 
always a non-local process. As the obtained solutions have sharp 
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